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ABSTRACT 

The non-Maxwellian K-distributions have been detected in the solar transition region and flares. These 
distributions are characterized by a high-energy tail and a near-Maxwellian core and are known to 
have significant impact on the resulting optically thin spectra arising from collisionally dominated 
astrophysical plasmas. We developed the KAPPA packag4E for synthesis of such line and continuum 
spectra. The package is based on the freely available CHIANTI database and software, and can be used 
in a similar manner. Ionization and recombination rates together with the ionization equilibria are 
provided for a range of k values. Distribution-averaged collision strengths for excitation are obtained 
by an approximate method for all transitions in all ions available within CHIANTI. The validity 
of this approximate method is tested by comparison with direct calculations. Typical precisions of 
better than 5% are found, with all cases being within 10%. Tools for calculation of synthetic line and 
continuum intensities are provided and described. Examples of the synthetic spectra and S'DO/AIA 
responses to emission for the K-distributions are given. 

Keywords: Techniques: spectroscopy - Methods: numerical - Radiation mechanisms: non-thermal - 
Sun: UV radiation - Sun: X-rays, gamma rays - Stars: coronae 


1. INTRODUCTION 

In astrophysics, the emitted radiation is usually the 
only source of information about the physical condi¬ 
tions in the emitting medium. Physical properties of 
the emitting plasma are then derived by analysis and 
modeling of the observed spectra. For a long time, this 
has been done under the assumption of a local, equi¬ 
librium Maxwellian distribution. This is done even if 
the emitting medium is optically thin and therefore per¬ 
haps not dense enough for the equilibrium to be always 
ensured locally. Such assumption is at best difficult in 
dynamic situations with particle acceleration, as e.g., a 
high-energy tail is difficult to equilibrate collisionally, 
since the collision frequency scales inversely with , 
where E is the particle energy le.g.. lMever-Vern(^l2007l l. 
iScudder fc Karimabadil (j2013f l argue that, in the case of 
stellar coronae, the assumption of the Maxwellian distri¬ 
bution should always be violated at heights above 1.05 
of the stellar radius. If long-range interactions are in¬ 
duced, e.g., by reconnection, wave-particle interaction, 
or shocks, the particles in the system can become corre- 
lated and do not have a Maxwellia n distribution (e.g. , 


Collierl 120041: IVocks & Mannl 120031 IVocks et alJ 

2008; 

Drake et al.l 20061: Livadiotis & McComasI 20091 

2010, 

20I3I: Pierrard & Lazaif 120101 IGontikakis et al.l 

2013: 

La,ruing et af.| |201 3f). Bather, the distribution exhibits 

a high-energy power-law tail. The K-distributions 

are a 


class of particle distributions having a near-Maxwellian 
core and a high-energy power-law tail, both of which are 
described by an analyti c expression (|Vasvliunasl Il968t 
lOwocki fc Scudd^lI983L , see also Sect. 2). The k index 
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has b een shown to be an independent thermodynam ic 
index (|Livadiotis &: McComasl[200l . l20I0Ll20IlL[20l^ in 
the g eneralized Tsallis statistical mechanics ('e.g.. lTsallisl 
rfm l2?M fLeiibTieH[2n?il [2n?il . 

The K-distributions can be derived analytically in case 
of a turbulent velocity diffusion coefficient inversely pro¬ 
portional to velocity. This has been shown for the 
plasm a in a suprathermal radiation field (|Hasegawa et al.l 
Il985f). for electrons heated by lower hybrid waves 
(iLaming fc LeDrill2007[l and for solar flare plasmas where 
the distribution function arises as a consequence of 
balance between diffusive acceleration and collisions 
(iBian et al.ll2014ll . 

Indeed, in solar flares, the /t-distributions provide 
a good fit to some of the X-ray spectra of coro- 
nal sources observed during partially occu lted flares 
(|KasDarova fc Karlickvl[20?M lOka et al.ll201.1l l. although 
a second, near-M a xwellian distribution is al so present 
(|Oka et al.ll2013[l . iBattaglia fc Kontail (|20I3f) used the 
AIA and RHESSI observations of flares to derive the 
distribution function in the range of O.I to tens of keV. 
These authors shown that the distribution derived in the 
low-energy range from AIA does not match the high- 
energy tail observed by RHESSI. A possible cause of 
this mis-match is the assumption of Maxwellian distribu¬ 
tion in the calculation of AIA differential emission mea¬ 
sures (DEMs), which may compromise the analysis, es¬ 
pecially if the high-energy tail is present and observed by 
R HESSI. _ 

iDzifcakova et al.l (1201 If) have shown that the n- 
distributions can explain the Si III transition-region line 
intensities obse r ved b y the S'OHO/SUMER instrument 
(jWilhelm et al.l Il995fl. especially i n the a ctive region 
spectr a (see also iDel Zanna et al.l (j20I4[) L iTesta et "m 
(j20I4[) inferred presence of high-energy tails from the 
ana lysis of Si IV spectra obs erved by the IRIS spectrome¬ 
ter (iDe Pontieu et al.ll2014ll . The K-distributions are also 
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routinely detected in the sola r wind (e.g.. iCollier et al 


1996 : lMaksimovic et al.lll997al lbl: [Livadiotis fc McComa 
201C : lT7e Chat et alJ 201llk The high-energy tails at 


keV energi es can arise as a co nsequence of coronal 
nanoflares (jGontikakis et al1l2013fl that are also able to 
produce the “halo” in the so lar wind electron distri¬ 
bution (|Che fc GoldsteinI 12014) . Furthermore, a claim 
has been made that the K-distributions were detected 
also i n the spectra of plan e tary nebulae (iBinette et alJ 
iwa iNicholls et a,bl MTl [mil IDonita et al.l 1201 fill , 
although this has been challenged as a possible ef¬ 
fect of atomic data uncertainties (IStorev et al.l 120131 
[Storey fc Sochi![201411 . The kappa-distributions are also 
one of the possible explanations of the non-Maxwellian 
Ho; prohles detected i n the Tycho’s supernova remnant 
([Raymond et alJl2010tl . 

Although the K-distributions were detected in solar 
flares, transition region and solar wind, their presence in 
the solar corona is currently unknown despite numerous 
attempts at their diagnostics. A diagnost ics of the high- 
energy elec trons have been attempted bviFeldman et ap 
(|2007ti and iHannah et al] (|20ir)H . iFeldman et al.l ( 20071 ) 
investigated whether the He-like intensities observed by 
SUMER could correspond to a bi-Maxwellian distri¬ 
bution with the second Maxwellian having a temper¬ 
ature of 10 MK. These authors argued that no such 
second Maxwellian is necessary. However, this anal¬ 
ysis was limited to Maxwellian distribution and did 
not include the effect of a prop er high-energy power- 
law tail. IHannah et al.l (|2010f l used the X-ray off- 
limb observations of the quiet-Sun performed by the 
RHESSI instrument (|Lin et al.l l2002h to obtain upper- 
limits on the emission measures as a function of k. How¬ 
ever, for temperatures of several MK corresponding to 
the solar corona, these upper limits are large and in¬ 
crease with increasing k. Direct attempts at spectro¬ 
scopic diagnostic s using EUV l i ne in tensities observed 
by Hmode/EIS (iCulhane et al.l 200711 were performed 
by iDzifcakova fc Kulinov^ 2010 ) and iMackoviak et al.[ 
(|2013f) . Indications of non-Maxwellian distributions were 
found using the O IV-0 V and S X-S XI lines. How¬ 
ever, such analysis was problematic due to large pho¬ 
ton noise uncertainties affecting weak lines, atomic data 
uncertainties and the possible presence of multi-thermal 
effects that would complicate the analysis. Therefore, 
even diagnostics using only strong lines will have be sup¬ 
plemented by a DEM analysis under the assumption of 
a K-distribution. Under the constraints of the current 
EUV instrumentation, such DEM analysis typically in- 
volv es many different elements and ionization st ages (see, 
e.g., iWarren et al.ll201^ IMackoviak et al.l[2014D . 

All this leads to a requirement of reliable calculation 
of synthetic spectra involving many different elements 
and ionization stages. In this paper, we describe the 
KAPPA package for calculation of optically thin astro- 
physical spectra that arise due to collisional excitation by 
electrons with a K-distribution. This package, allowing 
for fast calculation of line and continuum spectra for k- 
distributions, is based on the freely available CHIANTI 
atomic database and software, cur rently in version 7.1 
(jPere et al.lll997l : iLandi et al.ll2013ll . The manuscript is 
organized as follows. The K-distributions are described 
in Sect. |2j Synthesis of line spectra and continua are 
described in Sect. [3] and Sect. [H respectively. Section 
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Figure 1. The Ac-distributions with k = 2, 3, 5, 10, 25 and the 
Maxwellian distribution plotted for log(T/K) = 6.20 (top). Col¬ 
ors and linestyles denote the different values of k. Approxi¬ 
mations of the AC = 3 distribution in the low-energy range with 
a Maxwellian distribution according to I Livadiotis &: McComad 
ll2009h and lOka et"^ lf20il'l are shown in the middle and bottom 
panels, respectively. 

A color version of this image is available in the online journal. 

[5] describes the database and the software implementa¬ 
tion. Examples of the synthetic spectra and the AIA 
filter responses fo K-distributions are provided in Sect. 
[6l Summary is given in Sect. 0 


2. THE NON-MAXWELLIAN k-DISTRIBUTIONS 
2.1. Definition and Basic Properties 

The K-distribution of electron energies (Fig . 
[J) is defined as (e.g., lOwocki fc Scudd^ 119831 : 
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iLivadiotis fc McConiasl[200^ 


UE)dE = A, 






1 + 


E 


re+1 ’ 


( 1 ) 


{K—3/2)kBT 


where the =r(K + l)/(r(/t — 1/2)(k — 3/2)^'^^) is the 
normalization constant and fee = 1-38 xlO“^® ergs“^ 
is the Boltzmann constant. The ^-distribution has 
two parameters, Tg( 0,+oo) and k £ (3/2,-|-oo). The 
Maxwellian distribution at a given T corresponds to 
K —>■ oo. The departure from the Maxwellian distribution 
increases with decreasing n, with the maximum depar¬ 
ture occurring for k— :>3/2. 

While the most probable energy A^ax = (« — 
3/2)kBT/K is a decreasing function of k, the mean 
energy (E) = /2 of a ^-distribution does 

not depend on k and is only a function of T. 
Because of this, the parameter T has the same 
physical meaning for the K-distributions as the 
(kinetic) temperat ure for the Maxwellian dist ribu¬ 
tion. Additionally. ILivadiotis fc McComasI (j2009H and 
ILivadiotis fc McComasI ()201(1f l show that the T also cor¬ 
responds to the definition of physical temperature in the 
fra mework of the gene ralized Tsallis statistical mechan¬ 
ics (|Tsallislll988L l2009tl , and permits the generalization of 
the zero-th law of thermodynamics. Note that this fact 
permits e.g. the definition of electron kinetic pressure 
p = nekBT in the usual manner. 

Note also that the K-distribution is not the only pos¬ 
sible representation of a non- Maxwellian distribution 
with a high-energy t ail (e.g., iDzifcakova et al.l 120111 : 
iChe fc Goldsteinll2014f) . Nevertheless, its analytical ex¬ 
pression and a single additional parameter k make it a 
useful special case of an equilib rium particle distribu¬ 
tion associated with turbulence (iHasegawa et al.l 119851 : 
iLaming &: Lepril[2007l : iBian et al.ll2014ll . offering a rather 

straightforward evaluation of various rate coefficients as¬ 
sociated with radiative processes (Sects. [3]and|4]). 


2.2. Approximation by Maxwellian Core and a 
Power-Law Tail 

It is straightforward to see from Eq. © that in the 
high-energy limit, the ^-distribution approaches a power- 
law with the index o f — (k -|- 1 / 2). On the other hand, 
IMever-Vernet et © (| 199511 and ILivadiotis fc McComasI 
( 200911 showed that. in the low-energy limit, the k- 
distribution behaves as a Maxwellian with 


Tyi 


k-3/2 

- —T 

«: -I- 1 


( 2 ) 


The low-energy end of a K-distribution can indeed be 
well approximated by a Maxwellian, if this Maxwellian 
is scaled by a constant 


cm(«;) = A ^ 



(3) 


so that the two distributions match at the most probable 
ener p = (/t - 3/ 2)fcRT/K (Fig. [H middle] see also 
e.g.. iDziKhkov^ (|200^ . Fig. 1 therein). 


lOka et aP (1201311 attempted to approximate the core 
of a K-distribution with a Maxwellian at temperature Tc 


Tc 


^-3/2 ^ 

K 


(4) 


Such Maxwellian core has to be adjuste d (Fig. [1] bottom) 
by a scaling constant (|Oka et al.ll2013l Eq. (3) therein) 


c(k) = exp(l) 


T{k + 1) 
r(K- 1/2) 


.-3/2 





(5) 


so that the two distributions match at E = k^Tc. This 
approximation by a Maxwellian core leads to a worse 
match at very low energies A —>■ 0 (Fig. [1] bottom). 

These approximations suggest that the K-distribution 
can be thought of as a Maxwellian core at a lower tem¬ 
perature with a power-law tail. 


3. LINE INTENSITIES FOR THE k-DISTRIBUTIONS 

In the optically thin solar and stellar coronae, as well as 
the associated transition regions and flares, spectral lines 
arise as a consequence of particle collisions exciting the 
ions in the highly ionized plasma. The total emissivity 
Eji of a spectral line Xji corresponding to a transition 
/ —> j > i, in a A:-times io nized ion of the element X 
is usually expressed as (e .g., IMason fc Monsignori Fossil 
11994 [Phillips et aI1l200^ 


e,, = ^Apn{X+’‘) 


he Aj^ n(A+'=) 

\ji Ue n{X+^) n(X) 


AxneHH 


= AxGx,ji(T, ne,K)nenii , 


( 6 ) 


where hRi6.62xlO ergs is the Planck constant, 
c = 3 xl0^°cms“^ is the speed of light, Aji the Ein¬ 
stein coefficient for spontaneous emission, and n(A/'"^) 
the density of the ion -\-k with electron on the excited 
upper level j. In Eq. [SI the latest quantity is usually 
expanded in terms of the ionization fraction n{X'^^)/nx 
(Sect. 13.11) and the excitation fraction n(A/*"^)/n(A+'') 

(Sect. 13.2p . There, n(A+^) denotes the total density 
of the ion -\-k, and n{X) = nx corresponds to the total 
density of element X whose abundance is Ax, with nn 
being the hydrogen density. The function Gx,ji{T, Re, k) 
is the contribution function for the line \ji. The inten¬ 
sity Iji of the spectral line is then given by the emissivity 
integral of emissivity along a path I corresponding to the 
line of sight 


= J AxGx,ji{T, ne,K)nenHdl, 


(7) 


where EM = f rieRndZ is the emission measure of the 
emitting plasma. 

The CHIANTI atomic database provides the observed 
wavelengths Xji and the corresponding Einstein coeffi¬ 
cients Aji , while the electron density Ue is a free param¬ 
eter. To complete the synthesis of line intensities for the 
K-distributions, the relative ion abundance n{X'^^)/nx 
and the relative level population n{X'^^)/n{X'^^) must 
be calculated. This is detailed in the remainder of this 
section. 
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Fe 



Figure 2. Example of the ionization equilibrium for iron and the K-distributions. Relative ion abundances for Fe X-Fe XVIII are plotted 
as a function of T and ft and compared to the relative ion abundances for the Maxwellian distribution from CHIANTI v7.1 (full lines). 

A color version of this image is available in the online journal. 


3.1. Ionization Equilibrium 

A common assumption in calculation of the relative 
ion abundance n{X'^^)/nx is that of ionization equi¬ 
librium, i.e., that the relative ion abundance is not a 
function of time. Then, the relative ion abundance is 
given by the equilibrium between the ionization and re¬ 
combination rates. In coronal conditions, the dominat¬ 
ing ionization proce sses are the direct io nization and the 
autoionization (e.g.- lPhillios et ani2008[) . while the dom¬ 
inant recombination processes are radiative and dielec- 
tronic recombination. Since these processes involve free 
electrons, all of these rates depend on T and k (e.g. , 
Dzifcakov^ll992l:IAnderson et ai.lll996HDzifcakoy^l200^ 

Wannawichian et al.l 120031 iDzifSakova k, DudikI I2013D . 
We note that in the non-equilibrium ionization condi¬ 
tions, the n{X'^^)/nx depends on the specific evolution 
of the system, in particular o n the energy sources, sinks , 
and the resulting flows (e.g.. iBradshaw fc MasonI 1200.11 
iBradshaw et al.l l2004 lBradshawll2009D . Since radiation 
is an energy sink, the systern is then coupled. 

iDzifcakova fc DudiS (|2013D provide the latest available 
ionization equilibria for K-distributions for all ions of the 
elements with Z ^30, i.e., H to Zn. These calculations 
use the same atomic data for ionization and recombi¬ 
nation as the ionization equilibrium for the Maxwellian 
distribution available in the CHIANTI datab ase, v7.1 
l|Landi et al.ll2013l : lDerell200^ iDere et ^119971) . 

Figure [2] shows examples of the behaviour of the rel¬ 
ative ion abundances of Fe X-Fe XVIII with k. The 
ionization peaks are in general wider for lower k. Com¬ 
pared to the Maxwellian distribution, ionization peaks of 
the transition-region ions are in general shifted to lower 
log(T/K), while the coronal ions ar e generally shifted to 
highe r T, especially for low k = 2-3 (jDzifcakova fc DudikI 
[2f)ll) . Exceptions from these rules of thumb occur. E.g., 
the ionization peak of Fe XVII is shifted to lower T for 
K = 5, while for k = 2, it is shifted to higher T compared 
to the Maxwellian distribution. The shifts of the ioniza¬ 
tion peaks are typically Alog(T/K) Ri 0.10-0.15 for k = 2, 
although much larger shifts can also occur, e.g., for Fe VII 


(jDzifcakova fc Dudikll201,3l) . 

This behaviour of the individual ionization peaks with 
K strongly influence the resulting line intensities (Eqs. 
[6] and III). Therefore, the approximate temperatures de¬ 
termined from the observed lines in the spectrum can 
be different for a K-distribution and the Maxwellian dis¬ 
tribution. For a plasma where the high-energy tail or 
electron beams can be expected, the T is rela ted to the 
mean energy of the distribution (Sect. 12.11) including 
the high-energy tail. Notably, T can be very different 
from the Maxwellian “bulk” temperature Tm or Tq (Sect. 
E2l). Strong changes in the ionization equilibrium for 
the ^-distributions mainly in the transition region re¬ 
sult e.g. in the O IV being formed at log(T/Ar) r; 5.15 
for the Maxwellian distribution, but at R 5.0 for /t = 5 
and r 4.8 for it = 2 (jPudik et al.ll2014ah . The core of 
the distribution can have even lower temperatures - for 
\og{T/K) R 5.15, the log(Tc/K) = 5.0 for k = 5, but onl y 
4.55 for k = 2. The Tm are even lower (Sect. 12.2j) . 
Therefore, without a diagnostics of k in situations where 
the high-energy tail can exist, one has to be very care¬ 
ful in the estimation of the plasma temperature from 
the fact that a particular line is observed. The situa¬ 
tion is furthermore complicated by the dependence of 
line emission on the di fferential emission measure of the 
emitt ing plasma (e.g.. iWarren et al.ll2012l: iTeriaca et~^ 


20^), which is itself a function of k (jMackoviak et al 


20141) . 


3.2. Excitation Equilibrium and Rates 

The relative level populations n(A+^)/n(A+^) can be 
obtained under an a ssumption of excitation equilibrium 
(iPhillips et al.l[20081 Eqs. (4.24) and (4.25) therein). In 
equilibrium, the total number of transitions to and from 
any given level j is balanced by transitions both from all 
other levels m to the level j, as well as from the level j 
to any other level m. In the conditions of the solar and 
stellar coronae, ion-electron collisions are the dominant 
excitation mechanism, while deexcitations are facilitated 
either by spontaneous radiative decay (with the rates 
Ajm) and/or collisional deexcitation during ion-electron 
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collisions. The rates of electron excitation and deexcita¬ 
tion, and can be expressed as (|Brvansl [20061 : 
iDudik et al.ll2014^ 


/ 71- \ 
y/m^uJi \kBT J 


1/2 






y/^ujj V^bJ / 


( 8 ) 

(9) 


where ao = 5.29x10“® cm is the Bohr ra¬ 
dius, TOe = 9.1 X10“®® is the electron rest mass, 
Ih ~ 13.6 eV = 1 Ryd is the hydrogen ionization energy, 
oji and u}j are the statistical weights of the levels i and 
j, respectively, AEij = Ei — Ej is the energy of the 
transition, and Ei and Ej are the incident and final 
electron energies. The Tij{T,K) and Eji(T,K) denote 
the distribution-averaged collision strengths, given by 


Tii =Af^ --e 

k^T 


+ 00 


nj,{Ei) 


dE, 


^Ei- 


j^- -A 

" k^T 


+ C30 


+ (K-3/2)fcBT) 

) _ dAj 


^=+1 ae; 


^0) 


_ ^3 \ AEij 

/ (fc-3/2)feBT i 


( 11 ) 


In these expressions, ilji{Ej) = ilij{Ei) is the collision 
strength, i.e., the non-dimensionalised cross-section 




Ih tto ? 


E^ ot,{E,) 
‘ Ih 


( 12 ) 


where the a^i and erf are the electron impact excitation 


and deexcitation cross-sections, respectively. Note 
that with K^oo, the Tij{T,K) and Lji{T,K) revert 
to the Tii(T ] commonly us ed for the Maxwellian 


distribution (iSeatonI 119531 

Bnrffess & Tullvl 119921 

Mason & Monsienori Fossilll99'l 

: IBradshaw & Ravmondl 

2013|1. with the property of 

Ty(T) = Xji(T) being 


recovered. The Tij(T,K) and Eji{T,K), together with 
the Eqs. dS]) and m and the equations of statistical 
equil ibrium (Eqs. (4.24) and (4.25) in iPhillios et al.l 
[2?¥)1 . can then be used to synthesize the spectra 
for the K-distributions in the same manner as for the 
Maxwellian distribution (Sect. 15.21) . 


3.3. Collision Strength Approximation 

The calculation of the collision strengths for excita¬ 
tion and deexcitation averaged over K-distributions for 
large number of transitions introduces a problem of ac¬ 
cessibility of atomic cross-setions Ilji{Ej). Only a few 
database contain these data, and typically only for a 
small number of transitions. The CHIANTI dat abase 
and software (|Dere et al.l 119971 : iLandi et al.l I2013D con¬ 
tains spline approximations to the Maxwellian-averaged 
collision strengths for the majority of the astronomically 
interesting ions of elements H to Zn. CHIANTI allows 
for computation and analysis of solar spectra and is an 
important tool of the diagnostics of the solar plasma un¬ 
der the assumption of a Maxwellian distribution. 

We used the CHIANTI database to calculate the ap¬ 
proximate cross-sections H and subsequently approxi¬ 
mate excitation and de-excitation rates Ty (T, k) and 


Xji (T, k) for the At-distributi ons. This approxi mate 
method was describ ed e.g. in iDzifca koval (l2006ar) and 
tested for Ee XV by iDziRdkova fc Masonl (|2008ll . Here, 
we use this method to obtain the Tij(T, k) and Eji(T, k) 
for all transitions in all the elements and ions available 
within CHIANTI. 

The approximation works as follows: A func- 
tional form for the approx imation of H is assumed 
(lAbramowitz fc Stegunlll965ll 

^max 

n=J2 + D ln(n), (13) 

n—0 


where Ck and D are coefficients and u = Ei/AEij. The 
advantage of this approximation is a simple analytical 
evaluation of its integral over the distribution function. 
This approximation was often used fo r expression of the 
collision strength e.g. bv iMewd (fl^ . 

The Tij for the Maxellian distribution can then be 
written as: 


-y- _ ^Ejj 

“ k^T 


AE,-, /•“ 


e '“B 




which after integration leads to 

/ ^max \ 

T,, = Co + ^ yCn£n{v) + D£^{y) e\ (15) 


. fc=i 


where y = AEij/kAT and £n{y) is an n-th order expo¬ 
nential integral. 

The behaviour of H in the high-energy limit and the 
corresponding behaviour of T^- provide following condi¬ 
tions for the coefficients C„ and D for the electric dipole 
transitions 


dujjfij 
AE,j ’ 


^max 

Ty(-> oo) = Cn 


n=0 


£l{u 


1 ), 


(16) 


while for the non electric dipole, non exchange transi¬ 
tions 


Z7 = 0, T,,(^0) =Co, 

n-max 

Ty (-)> oo) = C„ = Vt{u = 1) (17) 

n—0 

and finally for the exchange transitions 
Co =Ci = H = 0, 

/ OO 

nd{u), 

^max 

Ty (-)> oo) = yy C„ = H(u = 1). (18) 

n—0 

The low- and high-energy limits Ty (—>■ 0) and Ty (—>■ oo) 
can be found in the CHIANTI database. The coefficients 
C„ and D can be evaluated from the collisional strengths 
in CHIANTI, averaged over the Maxwellian distribution 
by the least square method. To achieve the higher preci¬ 
sion, we used approximations (Eq. IT^ up to n-max = 7. 

The approximate method described here is also used 
to calculate the distribution-averaged collision strengths 
for deexcitation Eji(T, k). 
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0 IV 2s^ 2p - 2s 2p^ T,,, 



0 IV 2s" 2p "P ,/2 - 2s 2p" V,,, 



0 IV 2s" 2p "P ,/2 - 2s 2p" 'P 3/2 



Figure 3. Left: Approximation of T from CHIANTI. Middle: Comparison of the approximation of Q (red line) with atomic data of 
ILiang et al.l II2012I 1 without (black) and with smoothin (blue). Right: Comparison of the Ts calculated using our approximation (full lines) 
with direct calculations for the Maxwellian distribution (triangles) for different k distributions with k= 25 (violet), 10 (blue), 7 (turquoise), 
5 (green), 3 (yellow), and 2 (red). 


Fe XI Transition 1 -3 


Fe XI Transition 2-39 


Fe XI Transition 2-246 



Fe XI Transition 1 -3 




Fe XI Transition 2-39 





Figure 4. {top) and their relative errors {below) for the Fe XI 3s^3p^ ^P2-3s^3p'^ ^Po (left), 3s^3p^ ^Pi - 3s^3p'^(^D) ^So {middle), 
and 3s^3p'^ ^Pi - 3p^ 3d ^Fs (right) transitions. Black lines show comparison of the CHIANTPs approximation with the direct calculations 
for the Maxwellian distribution. Colors show the comparison of our approximation to direct calculations for the /^-distribution with /t = 25 
(blue), 10 (turquoise), 5 (green), 3 (yellow), and 2 (red). 


3.4. Validity of the Approximate Method 

Figure [3] demonstrates the approximation of O and 
calculation of T„ for the O IV transition 2s^2p^Pi/2- 
2s 2p^ ^P 3/2 at 1401.IhA. The atomic data for this transi¬ 
tion are taken from iLiang et al.l ()2012t) . We find a typical 
precision in the approximation of CHIANTI T’s of a few 
percent. This is the case for the O IV 1401. IhA transition 
shown, for which we find a precision of 1-2%. However, 
for a small part of transitions the precision can be signif¬ 
icantly worse, up to approximately 15%. 

Fulfilling the conditions (fT6ll - (|T8)) for the coefficients 
guarantees correct behaviour of H for high and threshold 
energies. It is however difficult to compare data for all 


transitions of each of ion. Occasional errors in the ap¬ 
proximation of O (Eq. 1131) cannot be excluded at present. 
Their propagation to the calculated of are further 
minimized by adopting 

'Y'approx 

y _ YCHIANTI /-jn'i 

-L K Maxwell ’ V-^^/ 

^ Maxwell 

where is the final T(k,T) for the K-distributions, 
T^ianti jg CHIANTI for the 

Maxwellian distribution, and and are 

T’s calculated from our approximations of the cross sec¬ 
tions for the K-distributions and Maxwellian distribution, 
respectively. 
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Fe XVII Transition 2-3 



Fe XVII Tronsition 2-3 



Fe XVII Transition 2-43 


Fe XVII Transition 3-50 






Figure 5. Same as in Fig. g] but for the following transitions in Fe XVII: 2s22p®3s lP 2 - 2s22p®3s ^Pi (left), 2s22p®3s lP 2 - 2s22p®4p 
^D 2 (middle), and 2s^2p®3s ^Pi - 2s^2p®4p ^D 2 (right). 


First tests of the precision of the approximate 
method desribed in Sect . 13.31 were performed by 
iDzifcakova fc MasonI (|2008[ ). These authors used u-max = 
5 and tested the validity of the approximation of the 
cross-section for some of the Fe XV transitions. An 
overall precision better than 10% was found. The ap¬ 
proximation worked almost perfectly for the alowed tran¬ 
sitions. Worse results were found for the forbidden tran¬ 
sitions. It was also found that transitions with strong 
resonance contributions and a low ratio of the excitation 
energy to temperature can also be problematic. How¬ 
ever, all the Hs for all tran sitions were reproduced to an 
accuracy better than 15% (|Dzifcakova fc Masonll200^ . 

To supp lement this analysis, we used rimax = 7 
(Sect. I3.3|l and tested the approximate method on two 
ions, Fe XI and Fe XVII. We used t h e orig inal atomic 
cros s sections f r om [ Del Zanna et al.l (120101) for Fe XI 
and IDel Zaniial (|201lh for Fe XVII. These Maxwellian- 
averaged Ty(T) are implemented in the CHIANTI 
database, version 7.1 (iLandi et al.ll2013ll . Here, we com¬ 
pare our approximation based on these Maxwellian data 
in CHIANTI with the Tij{T, k) and k) calculated 

directly from the Hs using the method of iDudik et al.l 

(IToTTbll . 

Figures. H] and [S] show several examples of the compar¬ 
ison of the direct calculation (hereafter, DC) with the 
approximate method for Fe XI (Fig. |4]) and Fe XVII 
(Fig. [5]). The DC are denoted by squares and the ap¬ 
proximate T«; by the full lines. Left columns in these 
figures show typical worst cases for strong transitions. 
We see that the error of the approximation depends on 
K and T; it typically increases with decreasing k. The 
worst cases are however still within 10% even for the ex¬ 
treme value of K = 2 considered here. Typical cases are 
shown in the middle columns of Figs. 0] and [SJ Here, 
the approximations are valid to within a few per cent 


for all KS. Finally, typical approximations for the weak 
transitions are shown in the right columns of Figs. |4]and 
O We again find that the approximations are valid to 
within f«10% for all ks. 

Figure [5] contain scatterplots of the relative error 
Tk/Tk, DC —1 plotted for each k at the peak of the cor¬ 
responding relative ion abundance. These scatterplots 
contain 447 transitions in Fe XI and 1050 transitions in 
Fe XVII that we were able to unambiguously indentify 
both in both the CHIANTI database and the atomic data 
themselves. The plots in Fig. |6] confirm that the approx¬ 
imate Tk do not depart from the directly calculated one 
Tk,DC by more than 10%. Typically, the relative errors 
increase with decreasing k; smallest errors are found for 
the Maxwellian distribution. Strong transitions typically 
have higher accura cy than the weaker ones, in a greement 
with the results of IDzifcakova fc MasonI (|2008[ 1. 

Finally we note that the approximation of T(k,T) to 
within 10% is considered satisfactory given the uncer¬ 
tainties in the atomic data themselves, which are typi¬ 
cally of the same order of magnitude, and the uncertain¬ 
ties of the spline-fits of the Maxwellian T{T) contained 
in CHIANTI, which are typically < 5%. 


3.5. Dielectronic Satellite Lines 

The rate coefficient for the dielectronic excitation from 
level i to level j and for an arbitrar y electron distribu tion 
funtion f{E) can be expressed as (jSeelv etTdlllQSTt) 


^diel _ 


meAEj, 


1/2 


WnTTiegi 


fiAE,,)Aa, ( 20 ) 


where gj and gi are statistical wieghts of double excited 
state and lower level, respectively; A®' is the autoion¬ 
ization (Auger) rate. The transition occurs at discrete 
energy AEji which corresponds to the energy difference 
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Figure 6. The relative error of to hdc XVII {right) as a function of at t emperatures corresponding 

to the maximum of the ion abundance ^K,(Tmax)' points are for the CHIANTI approximation (Eq. I19II . Different colors stand for 

results for ac = 2 (red), 3 (orange), 5 (green), and 10 (blue). 

between energy of states j and m. For the Maxwellian distribution, this equation leads to the well-known ex- 




























KAPPA: Synthesis of spectra for k-distributions 


9 


pression fe.g.. iPhillips et ani2008l Eq. (4.19) therein) 


C diel 


2(27rme*BT)3/2 g. 
For the ^-distribution, we have 




^diel _ 




9j 




2(27rmefcBr)3/2 g, 
which leads to 


1-H 


AEi. 


(K — 1.5)kBT ^ 


K+1 


C diel _ /^t 

^ — ^1 


diel 

Maxw " 


A^e '“B'r 


1 + 


AE,-. 


(K — 1.5)kBT 


K+1 


( 21 ) 


( 22 ) 


(23) 


4. THE NON-MAXWELLIAN CONTINUUM 

The continuum for the non-Maxw ellian K-distribution s 
is treated here using the approach of iDudik et al.l (I2012I1 . 
Contributions from the free-free and free-bound continua 
are considered. The two-photon continuum is not con¬ 
sidered, as its emissivity for K-distributions is not known, 
and its contr ibution is usually weak for the Maxwe llian 
distribution (|Young et al.l l2003t IPhillips et al.l I2008D es¬ 
pecially at higher densities. Nevertheless, at least for the 
Maxwellian distribution and a limited wavelength range, 
the two-photon continuum may not be a negligible con¬ 
tribution to the total continuum. We plan to implement 
it in the future. 


4.1. The Free-Free Continuum 

The total emissivity of the free-free continuum aris- 
ing due to electron-ion bremsstrahlung is given by 
(IDiidfk et al.l[2(rnl[2(rill 


A 

eff( A, K, T) = -- UeRH X 

xY,Kx{n,T)Ax 


9s{y,w) 




K+1 


dy, (24) 


where Ax is the element abundance relative to hydrogen, 
w = hc/XkBT = £/kBT is the scaled photon energy, the 
gs is the free-free Gaunt factor, and the Kx{k,T) is 
a function of k and T through the dependence on the 
ionization balance n^jux (see Sect. 13.11) 


Kz{^,T) 


1 32?! e® /27rfcB ^2 

471 3 meC^ V 3me ^ nz ’ 

k 


(25) 


where k is the ionization degree. The units of eg are 
ergs cm“^s“^sr“^A“^. 

The bremsstrahlung spectrum is strongly dependen t 
on K mainly at short wavelengths (|Dudik et al.l I2012D , 
where the tail of the ^-distribution strongly enhances 
the bremsstrahlung emission. Near the wavelength where 
the £ff peaks for the Maxwellian distribution, the free- 
free emission drops with k. At lar ger wavelengths it is 
enhanced again (see Figs. 2 and 3 in lPudik et al.l (|2012ll b 


4.2. Free-Bound Continuum 


on an excited leve l j is for the ^-d istributed incident 
electrons given by (|Dudik et al.ll20121 


efb( A, k,T) = 


£■5 


471 V tt (me/CBT)^/^ 


UeRH X 


E 

k,X 


^k+1 ^ gj _bf 

- Ax —cr An 

nx 9o 


1-f 


e-ii 


(k— 3/2)fcBT' ) 


K+1 


, (26) 


where £ = hc/X = E -\- Ij is the photon energy, Ij is 
the ionization potential from the level j with statistical 
weight gj, and is the ionization cross-section from 
the level j. 

A conspicuous feature of the free-bound spectra for 
the K-distributions are the greatly enh anced ionization 
edges (see Fig. 5 in iDudik et al.lPM!) . Generally, this 
increase comes from Eq. (1261) through the increase of low- 
energy electrons in a ^-distribution with respect to the 
Maxwellian distribution at the same T. However, details 
also depe nd on the ionizatio n equilibrium together with 
T and k (iDudik et al.ll2012Tl . 

5. THE KAPPA PACKAGE 

The KAPPA packagcQ currently allows for calculation 
of the synthetic spectra for integer values of k = 2, 3, 4, 
5, 7, 10, 15, 25, and 33, for which the ionization equi¬ 
libria are tabulated. These values should cover the pa¬ 
rameter space with sufficient density. The database and 
software for the KAPPA package is based on the IDL 
version of the freely available CHIANTI da tabase and 
software^ (iDere et a,l.lll997l : iLandi et al.ll2013[l . The rou¬ 
tines and database of the KAPPA package are contained 
in a standalone folder. It cannot be contained within the 
CHIANTI itself in order to prevent its automatic removal 
by CHIANTI updates. The path to the folder can be set 
by an IDL system variable in the idLstartup.pro file 

defsysv, '\datajpth'^ 'path to package'. 

The KAPPA folder contains the modified CHIANTI 
routines for calculation of spectra for the K-distributions, 
with the “dataJe” subfolder having the same structure as 
the CHIANTI’S “dbase“ subfolder. The modified CHI¬ 
ANTI routines follow the original CHIANTI routines as 
closely as possible. Their names end with an extra ” _fc “ 
before the .pro extension. The calling parameters of 
these routines are kept the same, except that the first 
parameter is always the value of k. 

The subdirectories within the database contain datas 
for ionization and recombination rates (Sect. 15. 1. 21) to¬ 
gether with the tabulated Tij(T,K) and iLji{T.,K) files 
in the ASCII format. Previous versions o f the mod¬ 
ificatio n corresponding to CHIANTI v5.2 (IDzifcakoval 
I2006b[l contained the coefficients for the approximation 
of £1. Then, the calculations for the ^-distributions were 
aproximately ten times longer compared to the CHI¬ 
ANTI for the Maxwellian distribution. Therefore, we 
decided to pre-calculate T{k,T) for a grid temperatures 
and K. These pre-calculated values of T(k,T) are con¬ 
tained in files names according to the ion and the value 
of K with the extension .ups, e.g., C-5-k2.ups for C V 
and k = 2. IDL saveffies containing the Tij(T,K) and 
Aji (T, k) are also provided. 


The emissivity of the recombination processes resulting 
in fc-times ionized ions of element X with an electron 


^ Iittp://kappa.asu.cas.cz 
3 www.chiantidatabase.org 




















































10 


Dzifcakova et al. 


At present, the KAPPA package fully corresponds to 
the atomic data contained in the CHIANTI version 7.1. 
Similarly, the routines provided in the KAPPA package 
are based on CHIANTI 7.1 routines, with the exception 
of routines for free-free continuum (see Sect. 15.3.11) . 

5.1. Ionization Equilibrium 

5.1.1. Ionization Equilibrium Files 

The ionization equili brium .ioneq and simil a r files were 
originally provided by iDzifcakova fc PudiH (I2013H . A 
minor software bug in the calculation of radiative recom¬ 
bination rates for the /c-distributions was found and cor¬ 
rected. This problem affected the ionization equilibria at 
log(r/ K) < 5 with the error being much smaller than the 
effect of K-distributions on the ionization equilibrium. 

These .ioneq files are produced in the same format as 
the original chianti.ioneq file. Therefore, these can be 
read by the CHIANTI routine read_ioneq.pro directly. 
The names of these files are kappa-02.ioneq and similar, 
where the numbers give the integer value of k. For more 
detail s on the .ioneq file format, see IDzifcakova &: DudikI 
(|2013ll . Appendix A therein. 

5.1.2. Ionization and Recombination Rates 

In addition to the ionization equilibria, total ioniza¬ 
tion and recombination rates are provided for each ion 
and a range of temperatures. Here, the total ionization 
rate is a sum of the direct collisional ionization rate and 
the autoionization rate. Similarly, the total recombina¬ 
tion rate is given by the sum of the radiative recom- 
binat ion rate and the total die lectronic recombination 
rate (|Dzifcakova &: Dudikll2013l l. These rates are stored 
in the respective database folder for each ion, e.g., the 
dbase/c/c-5/c-5-k25.tionizr is the total ionization rate 
file for C V and k = 25. The file format is ASCII. The to¬ 
tal recombination rate file has the same name except the 
.trecombr extension. The routines read-rate-ioniz-k.pro 
and read-rate-recomb-k.pro are provided for reading these 
files. 

5.2. Tools for calculation of line spectra 

The KAPPA package provides several routines for cal¬ 
culation of line intensities. These are listed in Table [TJ 
As already mentioned, these routines are based on CHI¬ 
ANTI routines, version 7.1. They can be used in the 
same manner as the CHIANTI routines, with the excep¬ 
tion that the value of k is always the hrst parameter. 

The most important of these routines is the 
pop-solver-k.pro routine that calculates the relative level 
population based on the distribution-averaged collision 
strengths 'Tij(T,K) and X jj(T , k) calculated using the 
method described in Sect. 13.31 Other routines for cal¬ 
culating line intensities (Table [T]) rely on this routine. 
Examples of synthetic spectra calculated for k = 2 and 
their comparison to the Maxwellian spectra at the same 
T are given in Sect. 16.11 

We note here that the method for calculation of the 
collisional electron excitation and deexcitation rates de¬ 
scribed in Sect. 13.31 cannot be applied to the collisional 
excitation by protons due to unavailability of the proton 
excitation cross sections. The proton excitation is typi¬ 
cally negligible, but may be important for some transi¬ 
tions. In the synthesis of line spectra, the proton excita¬ 
tion rate for K-distribution is assumed to be the same as 


for the Maxwellian distribution at the same temperature. 
It is currently unknown if this assumption is justified. 
Because of this, we advocate caution in using such lines 
for e.g. diagnostics of k from observations. 

An interactive widget for calculating the synthetic 
spectra is provided in the kappa.pro routine, based on 
the CHIANTI’S chss.pro. The value of k is selected by 
the choice of the ionization equilibrium. Subsequently, 
the excitation and line intensities are calculated for the 
same value of k. All other functionality of the chss.pro 
routine is retained. 

5.3. Tools for calculation of continuum 
5.3.1. Free-free continuum 

The CHIANTI database relies on the approxima- 
tions to th e Max wellia n bremsstra hl ung c alculated by 
lltoh et ahl (1200011 and iSutherlaiidl (1199811 and incor¬ 
porated in the freefree.pro routine together with the 
itoh.pro and Sutherland.pro routines. This approach can¬ 
not be followed in the modified CHIANTI, since no fit¬ 
ting formulae exist for the free-free continuum for k- 
distributions. Instead, we provide two options to cal¬ 
culate the free-free continuum for ^-distributions: 

1. Direct integration using the Eqs. (1241) and 
(l25l) . This approach is implemented in the 
freefree-k-integral.pro routine and requires the 
data-k/continuum/gffew.dat file cont aining the 
scaled gs{y,w) values provided by ISutherlaiidl 
(|1998ll . These gs values are then de-scaled and 
numerically integrated. Since the scaling depends 
on the ionization energy (and thus on the ioniza¬ 
tion stage k and the proton number Z), it has to 
be cy ried out for each ion separately (jDndik et al.1 
120121 1. Therefore, the direct integration using the 
freefree-k-integral.pro is time-consuming and im¬ 
practical. We note that, in practice, restricting the 
integration to elements with relative abundance of 
Az ^ 10“® introduces a relative error smaller than 
10“"^ and speeds up the calculations by a factor of 
«2. Note that the Maxwellian-integrated gff{y,w) 
are a part of the CHIANTI database. 

2. To overcome the long calculation times, the free- 
free continuum has been pre-calculated as a func¬ 
tion of Z for 101 logaritmically spaced temper¬ 
atures spanning log(T/K) = (4,9) with a step of 
Alog(r/K) = 0.05, together with 29 logarithmi¬ 
cally spaced points in A = 0.1 A-3 x 10^A with a 
step of log(A/A) = 0.2. These calculations are con¬ 
tained in tbe data-k/continuum/ff-kappa-02.dat hie 
and analogous hies for other values of k, each 
hie for a single value of n. The hies are in the 
ASCII format. The routine freefree-k.pro reads 
these hies, folds and sums them over the abun¬ 
dances to produce a semi-hnal free-free continuum. 
The hnal free-free continuum is then calculated 
for the user-input ranges of T and A within in 
the ranges specihed above. This is achieved hrst 
by linearly interpolating in log(r/K) and then by 
spline-interpolating in log(A/A). In this way, an ac¬ 
curacy of few per cent is achieved in the A = 1A - 
2 X 10"^A range, with the calculation time of few sec¬ 
onds. We note that this method should not be used 
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Routine name 

kappa.pro 
ch_synthetic_k.pro 
descale_diel_k. pro 

emiss_calc_k.pro 

freebound_ion_k.pro 

freebound_k.pro 

freefree_k.pro 

freefree_k_integral.pro 

isothermal_k.pro 

make_kappa_spec_k.pro 

plot_populations_k.pro 

pop_solver_k.pro 

read_ff_k.pro 

read_rate_ioniz_k.pro 

read_rate_recomb_k. pro 

ups kappa interp.pro 


Table 1 

List of routines within the KAPPA package. 

Function 

interactive widget for calculation of synthetic spectra, based on ch_ss.pro 
calculates line intensities as a function of k, n© and T 

converts Tij(T,K.) and Xj j(T, K)from the scaled domain _ 

for dielectronic satellite lines and performs correction in Eq. 
calculates hc/X Ajin{X'f^) 

calculates the free-bound continuum arising from a single ion 
calculates the free-bound continuum 

free-free continuum interpolated from pre-calculated data 
calculates the free-free continuum directly 
calculates isothermal spectra as a function of A 
routine for calculating the synthetic spectra 
calculates and plots relative level populations 
calculates the relative level population 

reads the pre-calculated free-free continuum as a function of Z and T 
reads the total ionization and recombination rates 
reads the total ionization and recombination rates 
routine for interpolating the T^^ (T, k) and A ji(T,K.) 


for calculation of free-free continua below lA and 
2 X lO^A, where the spline interpolation results in 
errors of several x 10 % or more. 

5.3.2. Free-bound continuum 

Using Eq. (l26l) . the free-bound continuum is straight¬ 
forward to calculate. The freebound_k.pro and free- 
bound_ion_k.pro routines can be used in the same manner 
as the CHIANTI’S freebound.pro and freebound-ion.pro 
routines. The only change is that these routines require 
a value of k as an extra input. Ionization equilibrium 
files (Sect. I5.1|l are read together with the cross-sections. 
The speed of the calculation is the same as using the 
original CHIANTI. 

6. SYNTHETIC SPECTRA AND 5DO/AIA RESPONSES 

In this section, we provide some examples of the 
calculated spectra that are of interest to the physics 
of the solar corona. Note that the behaviour of 
individual lines observed by the Hinode/FAS spec¬ 
trometer (|Culhane et al.l 1200711 and the possible ob¬ 
servational diagnostics with and without the effect 
of the ionization equili b rium are described els e where 
(iDzifcakova fc Kulinoval 120101: iMackoviak et al.l 120131 
iDzifcakova k, DudikI 120131 iDudik et ah 120141)11 , as is 
the a pplication for DEM diagnostics (|Mackoviak et al.l 

[Ml . 


6.1. Synthetic EUV Spectra 

An example of the synthetic isothermal spectra cal¬ 
culated using the isothermaLk.pro for the Maxwellian 
and K = 2 are shown in Fig. [Tj These examples show 
synthetic line and continuum spectra within the AIA 
171A and lOSA channels. The spectra are calculated 
for the electron density of Ug = 10® cm“^ and tempera¬ 
tures of log(T/Ar) = 5.9 for AIA 17lA and 6.2 for AIA 
I 93 A channels, respectively. Note that the temperature 
is kept the same for both distributions shown. These 
temperatures correspond to the maximum of the relative 
ion abundance of Fe IX and Fe XII under the Maxwellian 
distribution, respectively (see Fig. |2|). The line intensi¬ 
ties for K = 2 are decreased by a factor of several com¬ 
pared to the Maxwellian distribution. This is mainly an 
effect of the ionization equilibrium, with the maximum 


of the relative ion abundance for k = 2 being shifted to 
higher log(T/A") (Fig. [2l bottom). Note that the in¬ 
tensities of the continuum is several orders of magnitude 
smaller than the line intensities; therefore, the continuum 
is not visible in the linear scale on Fig. [71 
From Fig. [7] we see that at \og{T/K) = 5.9, the AIA 
171A channel is dominated by Fe IX independently of 
the value of k. Contrary to that, the situation for the 
AIA I 93 A channel and log(r /K) = 6.2 is more complex. 
This filter is dominated by Fe XII transitions between 
energy levels I--30 at 192.394A, 1-29 at 193.509A and 
1-27 at I 95 .II 9 A (iDudfk et al.ll2014bL Table B.4). How¬ 
ever, contributions from the 1-38 and 1-37 transitions 
in Fe XI at 188.216A and 188.299A are present as well. 
The relative contribution of these transitions to the total 
filter response to emission at log{T/K) = 6.2 increases 
from 4.4% for the Maxwellian distribution to 6.5% for 
k = 2. This is because this temperature is closer to the 
ionization peak of Fe XI than Fe XII for k = 2 (Fig. [7]). 


6.2. AIA responses for the k- distributions 


As an example of the usage of the synthetic line 
and continuum spectra we calculated th e responses of 
the Atmospheric Imaging Assembly (AIA iBoerner et al.l 
I 2 OI 2 I : iLemen et al.l l2012fl onboard the So lar Dynam¬ 
ics Observatory fSDO) IPesnell et al.ll20l'^ for the n- 
distributions. Note that even though the continuum 
intensities are weak compared to the line intensities 
(Sect. 16.111 at a particular wavelength, the continuum 
is a significant contr i butor to some o f the AIA bands 
(lO’Dwver et al.ll201(1l : iDel Zannall2013[l . This is because 
the filter response to plasma emission is given by the 
wavelength integral of the filter and instrument trans- 
missivity times th e emitted spectrum (e.g., Eq. (6) in 
IDudik et al.ll2009[l . 


The SDO/AIA responses calculated for K-distributions 
and log(ne/cm“^) = 9 are shown in Fig. |8l The peaks 
of the responses are typically flatter and wider, and 
can be shifted to higher log(T/K) for low k. This be¬ 
haviour is typical since it is given mainly by the ioniza¬ 
tion equilibrium (Sect. 13.11) . I t has been r e porte d for 
the TRACE filter responses b v iDudfk et al.l (l2009fl an d 
for Hinode /XRT responses by IDzifcakova et al.r ( 2012f) , 
where the AIA responses were also calculated for an 
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Maxwell, log(T/K) = 5.90 log(ng/cm ^) = 9 k = 2, log(T/K) = 5.90 log(nj,/cm ^) = 9 




Maxwell, log(T/K) = 6.20 log(n,,/cm ') = 9 k = 2, log(T/K) = 6.20 log(n,,/cm"^) = 9 



A [A] A [A] 


Figure 7. Example isothermal spectra at log(T/it') = 5.9 near the peak of the AIA 17lA wavelength response (top) and at log(T/K) = 6.2 
near the peak of the AIA 193A filter (bottom). The electron density assumed is log(ne/cm“®) = 9.0. 

A color version of this image is availatile in the online journal. 


earlie r set of atomic da ta corresponding to CHIANTI 
v5.2 (|Landi et al.l [200^1 . The AIA responses calcu- 
lated here repres e nt a s ignificant improvement over the 
iDzifcakova et ahl l)2012f) on es due to advan ces in the 
atomic data for AIA bands (|Del Zannall2013fl . 

We note that some of the secondary maxima, such as 
those at log(r/K) « 5.4 for AIA 17lA, 193A, 21lA, and 
335 A disappear for low k. This is again mostly be¬ 
cause of the wider ionization peaks and t heir relative 
contr ibutions to individual filter responses (|Dudik et al.1 
l2009r) . which gradually smooth out these secondary max¬ 
ima with decreasing k. We also note that the contribu- 
tion from Fe X a nd Fe XIV to the AIA 94A response 
(|Del Zannall2013l Fig. 3 therein) form a single, smooth 
secondary peak of the response for k ^ 5. 

These AIA responses for K-distributions can be used 
to obtain the DEMs using the regulariz e d DE M in- 
version developed bv iHannah fc Kontad (|2012H and 
iHannah fc Kont^ (|2013h . as well as for diagnostics of 
the distribution from combination of imaging and spec¬ 
troscopic observations (Dudik et al. 2014, in prepara¬ 
tion). 

7. SUMMARY 


We have developed tools for calculation of synthetic op¬ 
tically thin line and continuum spectra arising from col- 
lisionally dominated astrophysical plasmas characterized 
by a K-distribution. These tools constitute the KAPPA 
package, which is based on the freely available CHIANTI 
database and software. At present, the KAPPA pack¬ 
age can handle only values of k = 2, 3, 4, 5, 7, 10, 15, 
25, and 33, which should provide sufficient coverage for 
most spectroscopic purposes. Ionization and recombina¬ 
tion rates are provided together with ionization equilib¬ 
rium calculations. Approximations to the distribution- 
averaged collision strengths are provided. These are 
based on the reverse-engineered collision strengths ob¬ 
tained from the Maxwellian-averaged collision strengths 
available within CHIANTI. This is done for all tran¬ 
sitions in all ions available within CHIANTI, version 
7.1. We have tested the validity of this approximate 
method by comparison with the directly integrated col¬ 
lision strengths. For temperatures typical of the forma¬ 
tion of individual ions, typical errors of less than 5% 
were found. It was also found that the errors are always 
less than 10%. The errors are typically of the order of 
a few per cent for strong transitions, but the precision 
decreases for weaker transitions and low values of n. Con- 
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Figure 8. Responses of the AIA EUV filters for the K-distributions. Individual colors and linestyles stand for different values of k, as 
indicated. 

A color version of this image is available in the online journal. 


sidering the uncertainties in the atomic data calculations 
themselves, these errors are considered acceptable. 

Several routines for calculation of the synthetic line 
spectra, free-free and free-bound continua are provided. 
These routines are based on the CHIANTI routines and 
can be used in the same manner, except that the first 
input parameter is always the value of k. The calcula¬ 
tion of the free-free continuum is based on interpolation 
from pre-calculated values; however, an option of direct 
integration of the free-free gaunt factors is also provided. 
We aim to keep the database updated to reflect the newer 
releases of CHIANTI. 
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